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INTRODUCTION 


In an academic sense alone, a mathematical investigation 
of the stability of towed bodies and systems is of interest 
only if it adds to the existing knowledge in the general field 
of the course-keeping qualities of vessels. In a more practi- 
cal sense, it is an important and, in fact, necessary primary 
design aid which yields specific criteria to be satisfied for 
stable operation. These criteria are composed of terms which, 
for the most part, have distinct physical significance and 
which can be varied in design studies. Through an appropriate 
analysis of such terms, desirable operating characteristics may 
be found to a considerable degree of accuracy while such a de- 
Sign is still in the "drafting board" stage. 

This paper deals with the analysis of a system composed 
of a towboat or tug and a single barge. Both bodies are con- 
sidered to be unrestrained, hence the stability of the system 
will be dependent on the characteristics of the component 
parts. <A simple automatic control is specified for the towing 
body in order that the system should have directional rather 
than dynamic stability (or instability). 

In general the methods followed are those of Euler, Routh 
and Hurwitz and deal with systems of simple linear differential 
equations only. Laplace Transformations were not applied, 
since the author has found that many engineers are unfamiliar 


with the properties of this concept. 





MATHEMATICAL DEVELOPMENT 


The subsequent analysis presupposes the following basic 
simplifying assumptions: 

1. Both bodies have freedom in the horizontal plane 
alone, i.e., roll and pitch are ignored. 

2. The towline is assumed horizontal and taut at all 
times. 

3. Angular displacements are considered to be small, so 
that second order and higher terms may be neglected in the 
final equations of motion. 

With these premises in mind, consider a set of right 
hand orthogonal axes fixed relative to the earth and denoted 
by (x5-¥,)- Consider also similar systems of axes: (*,y) and 
(x,y) fixed in the barge and towboat respectively. Such a 
system is shown in Fig. 1 together with other symbols neces- 
sary to define completely the orientation of the system at any 


time. A complete and detailed nomenclature is given in Appen- 
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Kinematic Relations 
Certain expressions are obtained from the towing geometry. 


These are presented without comment based on Fig. 1 of Ref. 1. 


Yo, 92 V Sard — Ryde + x Sin) = me, 2 Sin (7 -W) 
x5, = J V Cos (p — Bde + x, Comp =>, - Cas (n - ¥) 
Vi = { V Sin (w — Bde + x Sim) (1) 


ip, ed V Cos ( — Ajde + x, Cos) 
.7=u+ are cos F Cos (Ww Pydt Jt Cos (W-Bjdt TF Cosw Eee Comb | 


yy 
=. + arc sin (¥ Sin (w — Fjde Al, Sin (U—Ajde + Sin Sin | 
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é-w= (2) 
P= T Gee Ty =—TSinf_ , Le ee Sint 
= = = 3 
i =—- | Cost Ty = T Sin” Looe x T Sin” (3) 


Differential Equations of Motion 


Since axes fixed in the bodies are essential to the de- 
velopment it is necessary to follow Euler, Ref. 2 in determining 
appropriate mass-acceleration terms. For the purpose of con- 
tinuity these expressions will now be derived. 

Choose a set of axes fixed in space and a set of axes 
fixed in the body as shown in Fig. 1. If P, and P, are the mo- 
mentum expressions relative to the axes fixed in the body then 
(according to Newton) the force components along the axes fixed 
in space are equal to the components of the time rate of change 


of momentum along these axes. Thus 





. = d ) - d ) /.. > . , 
y f - i a P. Cosy Py Sin | 
I eo | de 
= ( ) - « ) Se 4 ) ‘ 
%, a Poy de| P, Sanw Py Cost 
also 
X07 4 Gosw Y Sire! (53 
Y, ow Simy + Y Cosy 
._7= Pp. _ Py 
255 ) (6) 
Ys Py + Fe y) 
where 
Pre) u Py =M, uv : (7) 
Then 
YW, a -M, uw 
(8) 
Y=, vt My ou y) 


Y =e Tse ow 
; (9) 


Y=M,u+w, aw 
The moment of inertia terms are 

N=1U and W-Te (10) 
We consider here Ref. 1 and Ref. 2 and note that 


2M tk) | wt ky) 7 


where M is the absolute mass of the ship; /, is the moment of 
iments of the body about a vertical axis taken through its 
center of gravity, /, is the moment of inertia of the fluid 
displaced by the body taken about a vertical axis through its 


center of gravity and k,, k, and k are appropriate coeffi- 





cients of accession to inertia. Similar expressions may be 
written for the "bar" terms, thus 

Wo Wi(t+ ky) PSHE ek, ) | 

(az) 

We are now in a position to write the differential equa- 
tions of motion: 

l. Barge 

a) Longitudinal Force 


Ma —M, uw =X (u,v,W) + T, (13) 


b) Lateral Force 
M+ M, uwd=Y (u,v,h) +T, (1) 


c) Yawing Moment 


Iwo=W (u vv) + Noe) (15) 
ce Towboat 
a) Longitudinal Force 
Wa -W, 70 = 8 (20,05) + F, (28) 
b) Lateral Force 
8 +H, 00=F (0,0,0,3) +7, (17) 
c) Yawing Moment 
iv =N (a,%,¥,5) + Nery 18) 
3. The towboat rudder equations for rudder angle pro- 
portional to heading change are: 


a) Rudder order 


o =a, u 19) 


b) Constant time lag 
= 5 «eee um) 
dt 





Linear Simplifications 

Unlegr the previous SyStems are linearized, a classical 
solution i6 extremely difficult to obtain, in fact in a non- 
linear form the problem can best be treated through the use of 
automatic computing deviceu. However, only a closed form #olu- 
tion will be considered herein, hence the systems will now be 
linearized. 

T) to first agsumed that V ~ V , i.e. that both the tow- 
boat and barge are traveling at the same speed and that both ¢ 
and £ take on only small vaiues. Then the following relations 


are obtained from Figure 1 


i= VY Coal V 

vy = -V Sinfx~—-VF 

z= V Cosfex Y (21) 
v=-vV SinB x ye 


Noting these relations, the truth of the following iden- 


tities is obvious 


X (u,v,w) =X (Bw) 

y (u,v) = Y (6,W) (22) 
NO (u,v,w) = N (fh) 

C (a0,0.3) = ¥ (B¥,9) 

Y (a,0,0,3) = ¥ (Bu,3) 

N (a,0,0.3) = N (F,U,d) 


Limiting all displacements to the first order of small 
quantities, relations (22) may be expanded in Taylor's Series 


to obtain the following: 





<n 
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V(F,.0,d) =X, + ox (A — fy) + Mw -wy +H 7 -%,) 

of on 00 

os. Oe... ee: aw se 
F(B,U,d) = ¥, + 7 (A- BL) + =  - Uy + (8 - 2.) 

ee = ee a On ‘a — 
W(B,U,5) =, + (B-B) + ew - 0) +S 5-8, 


where the subscript, e denotes the equilibrium values about which 
the expansion takes place. Since it is assumed that initially 


the systemis on straight course andin lineitis obvious that: 


and 


y =X y= ¥,0+ Yuu 


(2u ) 


where the subscript notation has replaced the conventional par- 


Civgal derivative fowm. 





Continuing the linearization to include expressions (4). 


(2). and (3) the following relation# are established 
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ea ) far 
x x 
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In a similar fashion expressions (8), (9) and (10) become 
oe) 
Y =a Vb <i, OL 
i Set 
¥ =0 (26) 
Y=V,voe-Wve 
Nay 


With these linearizations the equations of motion become: 


l. Barge 


a) Longitudinal Force 


0 =a, + T (27) 
b) Lateral Force 
M Vu, V f= Vf + Vou (28) 


1 wa (u—L)dt -(t Pye] 
_T| =P y - <P 
7 : 





c) Yawing Moment r 


[U=NAf + Nyh — pT {v , (im Pde - fy (y —Fyde | 


~x T| +P <= & (29) 
Sh aw 7) 
°. Towboat 


a) Longitudinal Force 


i-85-7 (30) 
b) Lateral Force 


Mf vio -4H, y B= Yz6 Yau + ¥zo +T | wu ft] 
-1| ff (W—Bjde +7P n->P | (31) 
c) Yawing Moment 


Tp = Na+ Nl + Nz en T+ f¥ sya -- + (BF de | 


= 7) Py —2P 
e,7/ FP ia (32) 
d) Rudder Order 
8 =a, Wp (22) 
e) Constant Time Lag 
B= 5 + 7 
— (3h ) 


Equations (27) and (30) represent steady-state conditions 
and need not be considered further in the stability analysis. 
Non-Dimensionaiization of Equations of Motion 

Since the hydrodynamic forces may best be expressed in 


terms of dimensionless coefficients. one divides equation 


(25) by Sale ; Ml oe 
- O 9242 —— a4 PD 732 
(2) ays 7" : and (32) by S)'V 


Furthermore the following transformation is introduced for the 


/ 
independent variable ¢ $ a 


Here s i8 the number of tow lengths traveled, ¢ the true 
time, L the tow length and |} the linear velocity of the systen. 


Since V is considered substantially constant the following 


relation is written aot (35) 


9 





thus the following expre#sions are obtained: 
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UMing the furdgoing reiationg thy waqitationg of motion be- 


Barge 
a) Lateral Force 


a Pi mn PB -¥'far'gi' (37) 


‘ Ey f te - " (Pes B)} ees (fx +0N _- Li, |= N 


b) Yawing Moment 


‘ eof / Te 
ne Ub -N B= N wb (36) 
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Towboat 


a) Lateral Force 
~ ao 6G =e 
mn, £ Ww =§, #2 B - Yar - Yo 
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- ¥35 (39) 
- +; | fu + (E—£)} ds + oe es (132, i a 
b) Yawing Moment 
a’ f ?y - NB -Neh - N55 (40) 
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c) Rudder Order 
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d) Constant Time Lag 
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STABILITY ANALYSIS 


The equations of motion of the bodies under consideration 
together with the coupling effects due to the towline con- 
straint and the imposition of the control functions on the tow- 
ing body have been presented in the previous section. 

These equations are of general significance within the 
limits of linear theory with the exception that a constant time 
lag has been assumed for the control function. This assumption 
has been made purely to show that time lags of a sort can be 
included in the development. However in the subsequent analy- 
sis the control time lag s will be assumed to be zero. No 
loss in generality will occur since it is evident that the in- 
clusion of such a term will only serve to increase the order 
of the stability equation by one. Its exclusion will not make 
the analysis any less valid but will to some degree serve to 
Simplify it. 

To remove the integral signs in equations (37), (38), 
(39) and (lO). we operate on them by & and hence raise 
their order by one. The equations then become: 

l. Barge 

a) Lateral Force 


ie en!. Glee pe: i aaah i 
ai" vy | Flex “De "7k 
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(43) 
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b) Yawing Moment 


mee ard f / , e/ 
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a) Lateral Force 
of P] 
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~Ty +P p_¥38 = 0 
i y 
b) Yawing Moment 
’ SEE / —¢ J / ‘ 
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_ r) ; / / nays l / 4 
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(46) 
= apy a oT. me . 
- spy i R-N35 = 0 
? " 
c) Rudder 
S=aw (47) 


These are linear homogeneous equations in the dependent 
variables Y, oe W f B. 8, hence their solutions, excluding 
the highly improbable cases where the roots are equal, are 


given by the following equations: 


ce 


fas wer b=. 
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Be . ee ~ 48 
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It is our aim to determine: 

(1) the condition for which nontrivial solutions for the 
value 7, exist in the system 

(2) the condition that these 7, values are such that the 
motion be stable, i.e., that the real parts of all the %, values 
be less than zero. 

In order to fulfill condition (1) above, we substitute 
expressions (8) into equations (3) through (47) and realizing 
that the principle of superposition of solutions holds for 
linear differential equations, drop the summation notation and 
obtain the following relations: 

l. Barge 


a) Lateral Force 


b) Yawing Moment ; 
‘ raed ae fry / ‘ T 
jr 2 ie Nee. fo PS (@%,, Voy Bal 





a "4 
4 ol ,=t ae! / it 
u nN“ oO ioe B = *p pL ? (8 ; aps 4 uy 
B L pe L e t 74 l 
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2. Towboat 


a) Lateral Force 


- Whe 





b) Yawing Moment 
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c) Rudder 

oy a6 Wi (53) 


We now have a system of linear homogeneous equations in 
ds ie Vv, é, and 5. The condition that a nontrivial solu- 
tion exists among them is that the determinant of their coef- 


ficients be identically equal to zero. This means that 


ib ss see Ne Aq) sel 
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tion: 


wher: 
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Expanjion of this determinant leads to the following rela- 
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+ (aj + bh - Pf - md)(BU + EH + Le - MF - ND) 
+ (bj +eh-fg-mf- jg)(AJ+BH-LF - MD) (61) 
+ (ej amg = of Gai ly) 
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+ (ej - mg - jf) (Ed + Me - NF) 
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4 (g° f N+ g¢mM + gj L) (61) 
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+ (gomN + gg M) 
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To fulfill condition (2) on page (1,4), we need look to 
the theory of equations as developed for polynomials. In Hef. 
3, we find that the conditions for stability (the real parts 
of all the roots are less than zero) are stipulated by the 
following considerations: 

(1) all c, 70 (6: ) 


(2) all determintal relations obtained from 


aT Cy O 0 O 0 0 O 
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By striking out: 


first the first two rows and columns 


: > 
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second the first three rows and columns 
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i.e. CS Cy O 0 (67) 
C 3 Cy 
C C), C4 








ha 86 of Until. All Suth @e&préswiche have been written tr- 
cluding (6°). 

Traeea inégualities, when satisfied, indicate ttiability 
for th® System univer consideration. For any two plarticular 
bodi#s, a towline length and towpoint locations on both the 
barge and towboat may be tietermined to give desirable #tability 
@haratteristics. This, however, requires contid@rable numgeri- 
“ai ana.y#is as weli as time, unlews suitable computing fa- 
Rilities are available. 

The development of the theory as presented here is con- 
sidered to be a basis for future numerical investigation of 
the directional stability of any two bodies separated by a tow- 
line of arbitrary length attached to arbitrary points along 


the centerline of each body. 
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A SIMPLIFIED THEORY 


The foregoing analysis ha® been devoted to the direction- 
al stability of a tow@d configuration of two bodi®s, Separated 
by a towline of arbitrary length, attached at arbitrary points 
along the centerline of each of the bodies, and further con- 
fined by a control rule for the towing body proportional to 
its heading deviation. 

It. ig evident from this analysis that the results ob- 
tained defy, without considerable numerical calculation, ra- 
tional qualitative results. Thus in order to correlate quali- 
tative as well as some quantitative predictions with data from 
test measurements, certain rather broad simplifications will 
be made. 

It will be assumed that 

(1) the towline is horizontal, and taut (as before); 

(2) the barge does not roll or pitch 

(3) the towboat proceeds on a straight and prescribed 

course regardless of the motion of the barge. 

With this latter assumption, all the towboat functions 
vanish and only the barge functions remain. 

It is also convenient to define the transformation and 


dimensionless coefficients presented on pages 9 and 10 as: 
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With these in mind, equations (3) and (hh) reduce to: 


Lateral Force 


[vedi or [ eB] i op 


-m, B =e (8-76 = 0 (69) 
Yawing Moment 

cere A ee / xe of 
nb — N" oh 3 = yi [ 1 z 4] vy 


x’ T’ jo ee Te 
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Reasoning similar to that presented on pages 13 through 17 


leads 


where 


to the stability equation 

ec = 0 (71) 
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The motion will be stable if and only if all the C 


J 


values are positive and the following relations are satisfied: 
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In the case of the untowed body, it is noted (Ref. 1) 


that 
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2 (7, + %) 


where 7, and ©, are the dynamic stability indices of the un- 


towed body. 
Cy a my 
CS = -n'm, (o, 
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+ ao 
ae es / (1 7p) par Tip pn 
Tol =p = ”” a ?! Se Se 


/ ee / 
V'p-N'g | fT Pacer 


Then the oe values may be written as 


(76) 
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Without further mathematical manipulation, the following 


analysis may be made: 


be positive. 


1) Cy 


Then: 


>? 0 always 


by definition, all the C, values must 


4 
+ ns )—wN Ai 


~ 2 < 





2) 


a) 


b) 


c) 


oP 


a) 


b) 


¢) 


5 


> 0 Pom 


all dynamically stable untowed bodies, since 


and o ia) 8 


S 2 


l 
all dynamically neutrally stable untowed bodies, 
since 


c.. = O and o, <0 


all dynamically unstable untowed bodies for 
which 


Io,| < |lo,|,where 7, > O ando, <0. 


As a consequences of this, a towed body cannot 


be stable on course if |7,| > I7,| and 
IF,1 > 0. 
? 0 fom 


all dynamically stable untowed bodies, since 
the towing functions are positive and since 


both co, and o, are negative; 


all dynamically neutrally stable untowed bodies, 
since 


¢,=Oand G,< 0; 


all dynamically unstable bodies for which 


f 
x 4 
/ / y) ra |: / i — <a 

xy T my (1+ £ )+n'—t+n Moe, 0, ° 0 . 


> O for all towed bodies (see item 5 below), 


since by definition, (x) aa —N g) 7? 0, and usu- 


ally Grr =p, & He) > 0, ‘'y > 0 and T’and 


Fee 
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Cy > © by Gfinition. Since 


J al id - ri 4 / 
@,° T g= ®) i [Wao] -aa > 6, 


this condition may bé fulfilled @ither by increasing 

se (moving the towpoint farther forward), or. if 
this is in @ practical sense impossible, by increas- 
ing the drag in such a way that the static moment 
rate will remain unchanged or is decreased. The 
condition C. > 0 defines the location of the tow- 
point: it must always be a distance x’p forward of 
the center of gravity of the ship, which is greater 
than the distance between the center of gravity and 
the center of pressure of the static force Ya 

Assuming that the a values are all positive, it is ne- 


cessary to examine the determinantal relationships 
i= ane. bio to pl fal + ae: in’ 6. Gg 
a) , == AM, (7, +5.) Xp ( 7 } Mi, *n > mn oO, F, 
/ ‘ / fv / ee: Sey é / 7 / + f 
_n' mi [tp ¥’p-Na| ( BT's xg G-mi tay) 


Nit > 0 
2 


~an'nd (F,+3,) 5,3, - rye, +5) Q+ A | >O (77) 


where 





and 


i 


and Q, Rk, 


eee ai . de io -a/ 1A 
xp Va— Wall te] [ap emai mM» = wilt 


and n- n, must be positive. The following cases 


will now be examined: 


i) 


ii) 


Untowed body, dynamically stable. Here, 

~ Ose (co, + Gye, > 0 
always for 7, 0, > O and (0, + 7,) <Q. Thus, the 
sense of the inequality depends on the magnitude of 
T’R; hence, to achieve stability, it is desirable to 
keep this term small. This may be accomplished by 
proper selection of ie , the towpoint, and 7”, the 
towline length. 
Untowed body, neutrally stable. The first term 
vanishes since 0, 7, = 0; thus & must be made small 
in order to satisfy the inequality, and insure sta- 


bility. To accomplish this, the same reasoning is 


applied as in a). 


iii) Untowed body, dynamically unstable. Here, 7, > O 


and as long as G, < 0 and |@,| < |@,| then 

-n'm, (oc, “oyec, < 0 
Thus, for stability, it is desirable to make & as 
small as possible and to increase the drag in such a 


fashion that the inequality is satisfied. 
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If it is assumed that the conditions specified for 4, 
have been met, then the remaining factor to consider is the 
magnitude of Cc," (C) /C3). C, is determined by the values of 
Gite fand egy @ @nd 0 << O/C, << 1. It is therefore de- 
Sirable to make C, as small as possible, still keeping it 
positive, and thus insure the stability of the system. 

The criteria of directional stability will now be deter- 
mined for the case where the towline has increased in length 
to a point where the terms containing its reciprocal are 
negligible. Mathematically, this is achieved by letting the 
towline length increase without limit. In this case, it is 


founeetnat 


A, = =n'@, ras Yee v5 (x (7, + 7 Naked +n oil 
vl »/ / 


Since n, ns . T’ and (x, Y ‘, —WN -) are positive, this ex- 


pression will be positive if 
- = ‘ —— a7 a / / 
fm,’ (F, + F,)x5T +05, F) + (x5 ¥'g— Big) T <} O 


or 


(79) 





For all dynamically stuble ships, the right-hand side of the in- 
equality is positive; thus en‘ ig 4 —N‘,)T* must be kept as 
small as possible if the resulting motion is to be direction- 
ally stable. For dynamically unstable ships, the directional 
stability depends on e, T +n! 7, %,) > 0. Since in this 


case 7, 7,<0, it must be specified that for directional sta- 


i, Se =-_ = . is Bais 
balhity, Xp | ae cca 7, in addition to the condition that 


/ - = ; i 
5h, -N JT <—ap (6, +o,) (a, Tee, 2,) 


To investigate the consequences of decreasing the towline 
length to zero, i.e., either towing the ship froh its bow or 
with a bridle, we discard the terms containing the towline 
length » in equations (69) and (70) on page 22 since they have 
no significance in this situation. The a stability 


equation is then a cubic of the form of 


See 
> - Saag O (80) 
j=0 

where 
C= nat 
C, = (n Y p- ™ n’y) ™ ae ms (oa 2057) 
= / / =e 
Co, = (YY ym Ji, -) oh + my _ T = n'm, 7,9, +h, xp i 
Ca eg Mg) T (81) 


The criteria for stability specifies that all the C 
values must be positive and that the following determinant is 


satisfied: 





ar 


C) Co O 
A C t ral > : 

es C. C, O (62) 
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This reduces to (C, C,, Cy C,) O or 

anny? (F, +5) [no 3,424 1’|-n’my (2g ¥'g—W 4) 7 > 0 


Hence, the motion will be stable if 

=n,’ (o, + G,) pn" ao, +x, r’h (xp Yin —N - (83) 
For dynamically stable ships, the left-hand side of the inequal- 
ity is positive; consequently, the stability or instability of 
the towed ships is essentially a matter of towpoint location. 

For dynamically unstable ships, however, the extent of 
instability will be an important factor in determining the di- 
rectional stability, since the sign of the left-hand side of 
the inequality will depend on the sign of (no, 7, + ee 7). 
Thus, for the resulting motion to be stable, (no, 7, + ar T’)>0, 
Hence, the dynamic stability criterion 7, must be small if this 
condition is to be satisfied. 

Reducing towline length, however, is of very little prac- 
tical consequence, since an excessive decrease in towline length 
causes hydrodynamic interaction between the towing and towed 
Ships, a fact that will cause variation in the hydrodynamic 


"Constants". Thus, the basic assumptions of the analysis break 


down and the conclusions derived become untrustworthy. 
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NUMERICAL AMALYSIS 

The results of an investigation to determine the hydrody- 
namic coefficients of a typical barge form with various skeg 
configurations made on the rotating arm of the Experimental 
Towing Tank afford the data necessary for an analysis based on 
the foregoing treatment. 

The particular model considered had been found to be 
stable according to the free-towing standards set by the Experi- 
mental Towing Tank. These standards stipulate that if a barge 
is towed from its bow with a towline length of twice the barge 
length, that the barge should not be laterally displaced more 
than et of its beam. The towing device is a carriage that pro- 
ceeds along the tank rail at constant speed for the major portion 
of its travel. 

The accompanying sketch shows the bare hull of the barge 
model together with the skeg modifications investigated on the 
rotating arm (Ref. 6). The accompanying tabulation gives de- 
tails of the configurations as well as the comparable results 


of the standard free-towing stability criterion set by the Ex- 


perimental Towing Tank where these results are applicable. 


— a. 








one 


Model Particulars: 


Case I - 


Case II - 


Case III - 


Case IV - 


Case V - 


Case VI - 


00219 m. = 0.0351 n/ = 0.0219 
6.5 ft. x, = 0.508 2 = 2.0 
1.9 


bare hull (unsatisfactory by free-towing method) 


bare hull and skegs (5° to Z ) (unsatisfactory by 
free-towing method) 


bare hull, skegs and stabilizing prisms (satisfac- 
tory by free-towing method) 


bare hull and slotted skegs (not tested by free- 
towing method) 


bare hull and slotted skegs deflected to 6° (not 
tested by free-towing method) 


bare hull and slotted skegs deflected to 12° (not 
tested by free-towing method) 


In each case the model was tested on the rotating arm at 


Y 


a speed length ratio of a = 0.463 where Vis measured in 


knots. The following values of the hydrodynamic coefficients 


were determined: 












Case 


I | 0.008 
0.019 





0.0058 





| 









Calculation in accordance with the modified development 


(pg.21) and Ref. 1 leads to the following values for the untowed 


stability criteria (%,) together with the towed stability 





crit@ria 7, a& well as o pads, . 


CASE | case [| | CASE Lid 


+ 1 O47 wOleoe +. O°O735 
2. 188 2 131 


0 7687 | 0 76597 | 0 7°97 | 0.7697 
1 6314 





—0.0532 
—O O¥ae 


en oe 
0.044 || +0.4775i || +0.509i 
+0.248i | —O.1935 | —0.2405 | -—O 1735 
—1 873 | -—1 8155 | —1.8195 —1.9200 


eee ee ee 


The tabulation shows that none of the cases considered is 


0.82 





dynamically stable since in each case one of the stability 
roots (9,;) is positive. Of all the skeg configurations, only 
Case III produced a directionally stable barge when towed 
(since this is the only condition where all of the real parts 
of the values are negative). Case III was also considered 
stable as a result of the free-towing tests previously des- 
cribed. However, the drag of the barge in this case is notice- 
ably higher than for any of the others. 

It will be seen that none of the other skeg configura- 
tions creates such high resistance. While these other cases 
are not directionally stable, it appears that the instability 
may well be reduced without creating excessive drag. It does 


not seem too unreasonable to assume that further experim@utal 





investigation will yield tkeg forms that will produce direction- 
ally stable barges without adverse drag increments. 

The presences of complex 7, values (paired in the table) 
indicates that the motion is characterized by oscillatory modes. 
In all cases except Case III, the oscillation is undamped and 
tends to build up with time. This fact imposes continuous al- 
ternate yawing on the towed ship. In Case III, which is di- 
rectionally stable, the oscillations still appear but are 


damped and tend to diminish with time. 





CONCLUSIONS 


The preceding analysis shows that while the stability cri- 
teria and equation, for the case when two bodies are considered, 
have been determined, their evaluation for a given set of physi- 
cal and hydrodynamic constants is extremely laborious. It is 
difficult, if not impossible, to gain an extensive amount of 
qualitative information from the coefficients of the stability 
equation and the stability criteria because of the complicated 
way in which the various terms comprising these expressions ap- 
pear. 

In the simplified analysis where only the towed body is 
considered the following generalizations may be made. 

L. If a towed ship is to be directionally stable, the tow- 
point must be located forward of its center of gravity greater 
than the distance between the center of gravity and the center 
of pressure of the static lateral force. 

2. Depending upon its degree of stability, mass, towpoint, 
resistance, and moment of inertia, a dynamically stable ship may 
become unstable when towed with a long enough towline. 

3. Depending upon its towpoint, mass, moment of inertia, 
and damping functions, a dynamically unstable ship may become 
directionally stable when towed with a short enough towline, pro- 
viding that there is negligible hydrodynamic interaction between 
it and the towing ship. 


l). The motion of a towed body, which may be either di- 





rectionally @table or unstable, ip often characterized by ofcil- 
latory modes. 

5. In cases where a ship is extremely unstable dynami- 
callv, no posSible towing configuration can make it direction- 
ally stable. 

6. Increase in hull drag aft causes variation in each of 
the hydrodynamic parameters. Sufficient drag will cause most 
towed ships to become directionally stable. The experimental 
data presented indicate that the drag increments may be mini- 
mized and directional stability achieved by alteration and im- 


provement of skeg configurations. 





APPERMDIX ! 
NOMEN. LAT URE 


All lin@ar dimen#ion# are in feat an angular dimensions 


in radians. The time unit is thw #éecond. 


aw 
B= 
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Naas | 


m™ wy sl 


barge length 


barge towline point of attachment forward of barge 
center of gravity 


barge center of gravity 
towboat length 


towboat towline point of attachment aft of towboat 
ce@ater of gravity 


towboat center of gravity 
barge heading angle 

barge yaw angle 

towboat heading angle 
towboat yaw angle 

towboat rudder angle 
towboat rudder angle order 
tame 

towboat rudder time lag 
barge linear velocity 
towboat Linear velocity 
towlLine length 


over-all length of tow, i.e., from bow of towboat to 
stern of barge 


mass of barge 

mass of towboat 

longitudinal virtual mass of barge 
lateral virtual mass of barge 
longitudinal virtual mags of towboat 
lateral virtual mass of towboat 
momént of inertia of barge 

moment of inertia ofr towboat 


virtual mcm@nt of inertia of barge 





I virtual moment of inertia of toxboat 
i momént of inertia of diSplacedé fluid of barge 
» mom@nt of inertia of displaced fluid of towboat 


ar longitudinal virtual maws coefficient of barge 
k, lateral virtual mass coefficient of barge 
k, longitudinal virtual mass coefficient of towboat 
k, lateral virtual mass coefricient of towboat 
k’ virtual inertia coefficient of barge 
kb’ virtual inertia coefficient of towboat 
T towline tension 
re) density of water 
y hydrodynamic lateral force of barge 
N hydrodynamic lateral moment of barge 
Y hydrodynamic lateral force of towboat 
N hydrodynamic lateral moment of towboat 
’first, second, and third derivatives with respect 
to time 
’ indication of dimensionless coefficient 
subscript notation denotes a partial derivative 
with respect to the subscript; thus 
oy sy gee, CC. 
a: a? 
m, dimensionless longitudinal virtual mass coefficient 
of barge 
ms dimensionless laterai virtual mass coefficient of 
barge 
m, dimensionless longitudinal virtual mass coefficient 
of towboat 
m, dimensionless lateral virtual ma@s coefficient of 
towboat. 
n’ dimensionless virtual moment of inertia coefficient 
of barge 
n’ dimensionless virtual moment of inertia coefficient 


of towbout 
All other auxiliary svmbols are defined where they appear 
in the text. The exact expreszions for the dimensionless factors 


are given in each instance where they are first introduced. 


et 
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